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Abstract 



In this paper we construct the exact representation of the Ising partition function in 
' the form of the SLq{2, i?)-invariant functional integral for the lattice free (/, g)-fermion 

field theory {I = q = — !)• It is shown that the (Z, (7)-fermionization allows one to 
re-express the partition function of the eight-vertex model in external field through 
functional integral with four-fermion interaction. To construct these representations, 
we define a lattice (/, (7, s)-deformed Grassmann bispinor field and extend the Berezin 
integration rules to this field. At I = q = —1, s = 1 we obtain the lattice {I, g)-fermion 
field which allows us to fermionize the two-dimensional Ising model. We show that 
the Gaussian integral over (g, s)-Grassmann variables is expressed through the (g, s)- 
deformed Pfaffian which is equal to square root of the determinant of some matrix at 
q = ±l,s = ±1. 



^E-mail address: abugrijgluk.apc.org 



1 Introduction 



During the last decade, a considerable progress in understanding the meaning of the infinite- 
dimensional dynamical symmetries in the exactly solvable two-dimensional lattice spin mod- 
els and quantum field theories has been achieved. As it was shown by Belavin, Polyakov and 
Zamolodchikov in [1], conformal invariance of the critical fiuctuations in the two-dimensional 
lattice spin systems allows one to find a universality classes classification of the critical be- 
haviour of exactly solvable lattice systems and to assign to each class the corresponding 
two-dimensional quantum conformal field theory. This result is connected with the fact that, 
at critical points, such systems, besides the invariance with respect to the global conformal 
transformations {SL(2, C) group) possess larger symmetry, that is, the space of eigenvectors 
of the transfer matrix (Fock space for quantum conformal field theories) has an invariance 
with respect to the infinite-dimensional Virasoro algebra. 

However, as soon as the correlation radius of critical fiuctuations is finite nearby the crit- 
ical point, the conformal invariance is broken. Nevertheless, as it was shown by Zamolod- 
chikov in [2], exactly solvable two-dimensional quantum field theory can still possess in- 
finitely many integrals of motion, in the vicinity of the fixed point (this point corresponds to 
a quantum conformal field theory). This allowed him to suggest that the Fock space could be 
realized as an irreducible representation space for some infinite-dimensional algebra. This 
idea has been explicitly realized for exactly solvable massive quantum field theories. For 
example, in the case of quantum Sine-Gordon model, the Fock space has Uq{sl2) symme- 
try [3,4] and the deformation parameter q is connected with the coupling constant of the 
model. Also, as it was shown in [3,4] for the S'L'^(2)-invariant Thirring model, the creation 
and annihilation operators realize a representation of Uq{sl2) at g = —1. 

In [6] there have been also found infinitely many integrals of motion nearby the critical 
point of the eight- vertex model, corresponding to the two non-interacting Ising sub-lattices. 
In the remarkable series of papers by Jimbo, Miwa and their collaborators [7,8], it was 
shown that the eigenvectors of the corner transfer-matrix for the six-vertex model form the 
irreducible representation space of Uq{sl2) at — 1 < g < 0. Moreover, it was demonstrated in 
[9,10] that the space of eigenvectors of the corner transfer-matrix for the RSOS models and, 
in particular, for the two-dimensional Ising model [10] can possess this symmetry too. 

These results allow one to assume that the partition function of the two-dimensional 
exactly solvable lattice model can be expressed through the partition function of some 
5X5(2, C)-invariant lattice g-deformed field theory. 

In this paper, we make the first step examining this assumption and propose a sim- 
ple example of the construction of such a type (lattice field theory) in the case of the 
two-dimensional Ising model, making use of the exact representation of the Ising parti- 
tion function as the S'Lq(2, i?)-invariant functional integral in the lattice free (/, g)-fermion 
field theory {I = q = —1)- Moreover, it is shown that the method of g-femionization allows 
one to re-express the partition function of the eight-vertex model in external field through 
the functional integral with four-fermion interaction. 

For the construction of these representations we define a lattice (/, g, s)-deformed Grass- 
mann bispinor fieldQ and extend the Berezin integration rules to this field. Here / is a de- 

^Note that for / = e*'^^ the lattice (/, g)-deformed field one can consider as g-deformed non-abelian anyonic 
field 
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formation parameter for the commutation relations of "values" of this field in two arbitrary 
lattice sites, g is a deformation parameter for g-Grassmann spinor and s is a deformation 
parameter for commutation relations of two g-Grassmann spinors in a s-deformed bispinor. 
At/ = g = — l,s = lwe obtain the lattice (/, g)-deformed field which allows us to fermionize 
the two-dimensional Ising model (we call it the lattice (/, g)-fermion field). Moreover, we 
show that Gaussian integral over (g, s)-Grassmann variables is expressed through the (g, s)- 
deformed Pfaffian which at g = ±l,s = ±1 is equal to square root of the determinant of 
some matrix. 

It is a standard point that, in the scaling limit T — > Tc, a 0, the partition function of 
the two-dimensional Ising model is considered as the partition function of the free massive 
Majorana fermion theory [11]. One can obtain this result by representing the Ising partition 
function as the functional integral in the lattice real fermion field theory [12-16]. 

The fermionization method proposed in [13,15] proved to be effective for solving to the 
boundary condition problems and allows one to calculate the partition function of the two- 
dimensional Ising model with an arbitrary magnetic field attached to the boundaries [16]. 

Now we briefly recall the results obtained in [13,15]. Partition function of the Ising model 
on the torus is 



^exp 

{^} 



:i.r 



where r = {x,y) denotes the sites of the square lattice of size n x m, x = 1, . . . ,n, y = 
1, . . . ,m, /C = PJ, (jfl = ±1, and x, y are unit vectors along the horizontal X and vertical 
Y axes respectively. This partition function can be presented as the sum of the functional 
integrals over the lattice real fermion field 



Z = const / Vipexp 



where 
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Vxi^r = A+x, Vylpr = 1pr+y, t = tauh/C, Vlp = JJ'^V'r 



i,r 



[1.2) 

:i.3) 



:i.4) 



Here, for example, the indices P, A in Z^^ denote periodic and antiperiodic boundary condi- 
tions for the field ipr along the horizontal X and vertical Y axes respectively. The appearance 
of the four-component Grassmann (bispinor) field is related to the requirement of lin- 
earity of the quadratic form in ( |1.3| ) with respect to the shift operators Vx , V^. Note also 
that for such a quadratic form the doubling problem for the lattice fermions, arising in the 
course of usual formulations of lattice fermion field theories [17] is absent. 
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In the momentum representation, the matrix D has the following form 





-1 



-1 



1 



-1 

-1 - te-'P^ 



1 + te'P^ 1 
1 1 + te'P^ 







1 




(1.5) 



where for the periodic and antiperiodic boundary conditions along the X, Y axes, pi,P2 are 
integer and half-integer respectively. 

Using momentum representation, it is not hard to calculate functional integral ( |1.3| ) and, 
for example, one can obtain for the Z^^ 



Z^^ = consty det D = c JJ(s + s — cospi — C0SP2) 



:i.6) 



P1.P2 



:i.7) 



where s = sinh2/C and c = const, sinh2/Cc = 1 at T = T,,. 
Introducing the notation 

cosh 7(^2) = s + — C0SP2 
and calculating the product over pi, one gets 

Z'^^ = cY[ cosh -7(p) = c n exp(-ma^) (1 + exp(-ma^)) . (1.8) 
p 2 p 2 d 

Here a is a dimensional lattice spacing, and we define the dimensional momentum by sub- 



stitution P2 pa. In order to study the continuum limit, one expands ( |1.7|) into a power 
series in p around p = 0. Restricting ourselves to the quadratic terms in p, we obtain 



I'ip) 2(. - 1)^ 



1.9) 



Considering different limits in ( |1.9|) , one gets e{p) either for the massless or for the massive 
Major ana fermions: 

1). T = Te, s = l, a^O, 



/2s - 1 



fi = lim 

s = l V s a 



and e^(p) = p^ 



2). T ^Tr 



s ^ 1 , a ^ , 
= lim \ 



V s a 



\p2 and ^^{p) = + p^ ■ 



The free energy of the Ising model in the continuum limit turns into the free energy of 
Majorana fermions at finite temperature. To see this, let us define, a,sa-^{),n^oo.,m^ 
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oo, the system volume as L = na and the inverse temperature of the thermostate in which 
Majorana fermions are placed as /? = = ma. It is not hard to show that in the 
thermodynamic limit (L — > oo) the main contribution to the free energy comes from Z^"^. 
Then, for the free energy density / one obtains from ( |1.8|) 

oo 

/ = J rfpln[l+exp(-/3£(p))]. 

— oo 

This expression confirms (" +" in front of the exponential) that the partition function of the 
two-dimensional Ising model on the strip of width (3 = ma in the continuum limit describes 
the Majorana fermion system in the thermostat. 

This paper is organized as follows. In section 2 we give a short introduction to g-spinors 
and the quantum matrix group SLq{2, C). In section 3 we define the lattice (I, q, s)-deformed 
Grassmann bispinor field. In section 4 we define integral calculus for this field. In section 5 
we propose a method to represent the partition function of the two-dimensional Ising model 
with the nearest-, next-nearest-neighbour and four-spin interactions (eight- vertex model in 
external field) in the form of the functional integral over the lattice (/, g)-fermion field (/ = 
q = —1, s = 1). To realize this representation we use the GLg{2, C)-covariant generalization 
of the Berezin integration over Grassmann field proposed in section 4 and calculate the 
functional integral for the Ising model with the nearest-neighbour interaction. 



2 q-Spinors and the quantum group SLq{2, C) 

Let us briefly recall the known facts about g-spinors, their g-Grassmann realization and the 
quantum matrix group SLq{2, C) [18-23]. By g-spinor v we mean the two-component object 



(2.1) 



with the following commutation relation for the components 

o . 1 
-■i 

Q 

where g is generally a complex number. At g = 1 it is usual spinor. Following the paper 
[19], one can look at g-spinors as belonging to the quantum plane A^'". 

One can consider this quantum plane as the two-dimensional representation space for 
the quantum matrix group SLq{2,C) (g-deformed SL{2,C)) [18-23]. Using the i?-matrix 
for this group, one can rewrite commutation relation ( |2.1| ) in the covariant form 



(2.2) 



where Greek indices go over 1 , 2 and Rq is a symmetric 4x4 matrix (the rows and columns 
are numbered as (11), (12), (21), (22)) 



^11 p 



TDlP 







g-( 
1 




-1 





1 

g 



(2.3) 
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This matrix has eigenvalues q and —q ^ and satisfies the quantum Yang-Baxter equation 
[18-20] 

Rql2 Rq23 Rql2 = Rq23 Rql3 Rq23 , (2-4) 

and the Hecke relation 

Rl-{q-q-')Rq-i={i + qRq) (/ - q'^Rq) = 0, (2.5) 

where / = I^j^ = 5° 5jf is the unit matrix. One can write the following representation for 
the ^-matrix [20-22] 

Rg = qSq - q'^Aq, 

where matrices 



1 

1 + g 
1 

l + q 

are orthonormal projectors (to the subspaces with eigenvalues —q^^ and q respectively) 



^^ = VT^2(^-^~'^^)^ (2-6) 



Sq = y—^{i + qRq) (2.7) 



A ■ S = S ■ A = A? = A S"^ = S 

These projections are the quantum counterparts of the classical (g = 1) antisymmetrizer and 
symmetrizer for the tensors with two indeces. 

Using the g-antisymmetrizer Aq, one can rewrite commutation relations ( |2.2| ) 

The classical [q = 1) Grassmann spinor can be defined through 

<Jq ^p^ ^ '^i 

or, in components, 

Analogously, one can define the g-Grassmann spinor ip dX q^ ^ —1 [21,22] 

Sf^p r V = 0, = -q Rqi (2.8) 

or, in components, 

(^1)2 = (^2^)2 _ ^2 ^1 _ ^2_ ^2.9) 

It implies (see the paper [19]) that ip belongs to the quantum plane A^^"^, which we denote 
as ^(2). 

In order to show that commutation relation (|2.8| ) is invariant with respect to the trans- 
formation of by a quantum matrix A G SLq{2, C) 

r^^^^{%%)K)^{:)^K)^ (2.10, 
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we recall some properties of SLq{2, C). 

The matrix elements commute with the components if)"' and belong to the associative 
algebra of functions on the quantum group SLq{2), which we denote as A or FuUq^SL (2)) 
[20]. This algebra is a Hopf algebra. It implies that the following maps are defined in this 
algebra: 

a) comultiplication A 

where symbol " 0" denotes the tensor product of the quantum spaces, 

b) counit e 

where C is a complex number, 

c) antipode i 

A^A : ^(A^) = (-gf-^i^, 

where is the quantum minor of the matrix element A'^ defined for the quantum matrix 
A G SLq{2, C) in the following way [20] 

41-42 42_4i 4i_4i 42_42 

d) multiplication m 

A^A^A: m(A^® = A^A^, 

and the matrix elements A'^ obey the commutation relations 

^q^pA'J^A'^ = A'^A^R^j^j^. (2.11) 

which are equaivalent, due to ( p.lO|) , to 

ab = qba, ac = qca, bd = qdb, cd = qdc, be = cb, ad — da = {q — q^^)bc . 

These maps allows one to define the following operations with quantum the matrix A G 
5L,(2,C): 

a) A(y4) = A® A defines the rule of multiplication of quantum matrices, and matrix 
elements of A® A have the form A" (g) ( ® denotes the tensor product along with the 
usual matrix multiplication), 

b) e{A) = 1 defines the unit matrix in SLq{2,C), 1 = 6^, 

c) i{A) = A~^ defines the inverse matrix, that is i{A) ■ A = A ■ i{A) = 1, where "■" 
denotes the matrix multiplication, 

d) m{A ®A) = A^A defines the usual tensor product for quantum matrices, and matrix 
elements of A^A have the form A"^ AJ. 

e) the quantum determinant of the quantum matrix A 

(iet,i = 5:(-g)^-°A|i^ = l. 

a, 13 
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Using ( CT) , one immediately shows that it commutes with all the matrix elements A'^. 
Linear transformation ( |2.1CI| ) of the quantum plane "^^"^^ can be rewritten through the map 5 

^(2) ^^^^(2) . ^(^^^ ^ ^^^^ ^^^a) =A'^0 i)^. 

It means that \E'^^'' is a comodule of SLq(2,C). Then, using the map the action of 

quantum matrix A on the tensor product \E'*^^^ ® ^1'^^-' can be defined as the action of the 
tensor product A^A 

^(2) ^ ^(2)i®i(i^^(2)) ^ (i5^^(2)) , 

It is convenient to define the matrix A^A in the following way 

A^A = Ai-A2, A, = Am = A"pS;, A2 = mA = 6^AJ. (2.12) 
Hence, its action on \E'*^^) '^^'^^ can be written as 

{A^A)(^{iIj (^iP) = Ai- A2®{ij = 
{Am) ■ ® V) = (^®^) ® = ® (2.13) 

After using the matrices Ai and A2, commutation relation ( p.ll| ) can be represented in the 
form 

Rg-Ai-A2 = Ai-A2-Rg. (2.14) 

Let us show that quantum matrices in ( |2.10| ), ( |2.11| ) form the quantum matrix group 
SLg{2,C) with respect to the multiplication 0. It is easy to show that matrix elements of 
A" = A'^A satisfy commutation relation (|2.11|) , if A', A G SLq{2,C). Indeed, one defines 
matrices 

A[ = A2 = mA\ 

which satisfy the relations 

Rg ■ A[ ■ A'2 = A[ ■ A'2 ■ Rg , A'l = A[®Ai, A'^ = A'a^Ag. (2.15) 
Using the comultiplication property 

{A[ ® Ai) ■ {A^ ® A2) = iA[ ■ A^)(g){A^ ■ A2) 
and (121^) , (|2J[5| ), we get 

Rg ■ A'; ■ A'^ = Rg ■ {A[^Ai) ■ iA'^(bA2) = Rg ■ {A[ ■ A'^)®iAi ■ A2) = 

{A[ ■ A'^)^{Ai ■ A2) ■ Rg = {A[^Ai) ■ {A'^®A2) ■ Rg = A'l ■ A'^ ■ Rg. 

Hence, one concludes that the matrix elements of A" do really satisfy relation ( p.ll|) . More- 
over, the g-determinant has the property 

detgA" = detg{A'(§A) = detg{A') ® detg{A) . 

Then, the matrix elements A"f^ G Fung{SL (2)). 
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Now one can easily show the invariance of commutation relation (TB) under transforma- 
tion ( p.lOj ). To see this, it is convenient to rewrite ( p. 81) in the form 

(/ + qRg) . ® ^/>) = 0. 



Then, using (|2.13|) and ( |2.14| ), one proves the invariance by means of the following chain of 
equalities 

(/ + qRg) -(ip^ip) = (i + qRg) ■ (A^^tp) ® (A^ip) = 
(/ + qRg) ■ Ai ■ A2®iiJ ® ^) = ■ ^2 ■ (/ + qRq)®{ip ®i)) = 

{A®A)®{i + qRg)-{'4)®^) = Q. (2.16) 

3 Lattice (l,q,s)-Grassmann bispinor field 

In this section we define the lattice (l,q,s)-deformed Grassmann bispinor field which is needed 
for the (/, g)-fermionization of the two-dimensional Ising model. For definition of this field 
we use the results of the paper [24]. Recall that / is a deformation parameter of the com- 
mutation relations of this field at two different lattice sites, g is a deformation parameter for 
g-Grassmann spinor and s is a deformation parameter for the commutation relations of two 
g-Grassmann spinors in s-deformed bispinor. 

At first let us define (g, s)-Grassmann bispinor ip = {ipf} ( a = 1, 2 and i = 1,2) 

^-C:)- ^-CD- ^-Ci)- 

where ijj belongs to a quantum vector space and e \ ip2 e ^^2^. Then, from (2.8) 
it follows that ipi and ip2 satisfy the commutation relations 

C®^f = -9<4V^7®V'f. (3.2) 

We consider "^^^^ as the Hecke sum with respect to the matrix Qq [24] 

^(6) = ^(2)0^^(2) _ (3.3) 

This matrix defines commutation relations of ipi and ip2 in s-deformed bispinor ip. Following 
definition of the Hecke sum [25,26], let us consider the Q^-matrix as linear map 

Qg : — ^ ® ^ : 

Q,(^^ ® ^^) ^ Qf^^ ri ® ^2 = -^^2 ® V^f , (3.4) 

where in general case the deformation parameter s is a complex number. 

Following (3.2) and (3.4), we define the i?-matrix for \E'*^''^ as linear map 
R : (g) ^C*) > (g) ^W; 

^(C ® ^P^) ^ Rf^tl ® < = -X.t^; ® , (3.5) 
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where Xg is a diagonal matrix 



j-apkl 1 rfc d r ^ ra rfc d / r N 



and require that the i?-matrix satisfies the quantum Yang-Baxter equation 

R12 R23 R12 — -R 23 12 -R 23 , 

and the Hecke relation 

— (J's — is)R — Js — {Js -\- Js ' R) ' {Js — is ' R) — 0, 
where Jg — Js^pi\ — I^p Itj ^^i^ matrix, Jg is a diagonal matrix 

and Jg^Tg- Jg- 

It is not hard to check that the following anzats for i?-matrix 



3 ' 



where 



l,iik<l 
0/iik>l 



satisfies the Hecke equation (3.8). 

In the product '^^'^^ (8) ^^^^ in (3.5) we imply the following ordering: 

^(6) (g) ^(6) ^ ^(2) ^ ^(2)0^(2) ^ ^(2)0^(2) ^ ^(2)0^(2) ^ ^(2) 

With this ordering the i?-matrix has the form: 



R 



Rq 











[s-s-^)i Qg 



'9 







Rn 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



(3.10) 
(3.11) 



(3.12) 



(3.13) 



By analogy with (2.6) and (2.7), one can define from (3.8) the quantum (g, s)-anti- 
symmetrizer 

M^^'^ = {Jg-Ig-R), (3.14) 
and the quantum (g, s)-symmetrizer 



S^'^'^^^il + Jg-R) 



(3.15) 



which are orthogonal projectors 



^(q,s) . ^ {q,s) ^ ^ (q,s) . ^{q,s) ^ (^(9.^))2 ^ . (X, + X) ' 

The quantum symmetrizer 5*^'''*'-' allows one to define the commutation relations for the 
{q, s)-Grassmann bispinor {q^ ^ —1) 

5(^'^)-(V^®V') = 0, or ^^®iP^^+{X-R)tti^k®Vi (3-16) 

Let us require that this commutation relations are covariant with respect to the transforma- 
tion Tp hj & quantum matrix A e SL g{2, C). Then taking into account (2.4), it is not hard 
to find the following solution of the quantum Yang-Baxter equation (3.7) 

Qq^Rq. (3.17) 

Substituting this solution in (3.16), we get 

i^? iPf = -q iPf iPf , (V'f )' = 0, 

^2(g,^l ^ _s^i^^2^ (3.18) 
■02 ® V'l = -s il^l^i^l- s{q - q~'^)il^\ ® ill . 

To define the lattice (/, g, s)-Grassmann bispinor field ipf{r) we consider the d-dimensio- 
nal hypercubic lattice. The lattice sites are labeled by the vector r = (pi, p2i ■ ■ ■ , Pd), where 
Pi are integer numbers and the total number of sites is A^. For the sake of convenience, the 
lattice spacing a is fixed to be equal to unity. Let '0f(r') e "^f) be the (g, s)-Grassmann 
bispinor on the site r', which we consider as a "vahic" of the (I, q, s)-Grassmann field '4^f{r) 
on the site r'. The definition of the lattice (/, g, s)-Grassmann bispinor field requires some 
information about its commutation relation in two arbitrary sites r.^ = {li, I2, . . . , Id) and 
i^n = {ki, k2, . . . , kd). In this paper we suppose the quasi-one- dimensional ordering for the 
lattice sites: > r„, if Zi — /ci > and the others — ki are arbitrary; if li — ki — 0, then 
T^m > fn-i if ^2 — ^2 > and so on, and ri < r2 < < ... < rjv (m, n — 1, N). 

Let us denote a quantum vector space which is generated by 4A^-component (/, q, s)- 
Grassmann vector ip°'{r) as ^ {N is the number of sites on the lattice). By analogy with 
(3.3), we must consider the quantum vector space ^' as the Hecke sum with respect to some 
matrix Q 

^ = *We^*g)e^ . . . eQ^rW_^eQ^'S . (3.i9) 

For definition of this sum and its "^-matrix let us consider Q-matrix as linear map 
Q : * W ® *W *W * W 

»^ In ' m ' m ' n 

g(C(r„) ® i^ivm)) = QfJI ® iUrm) = -}c(^m) ® V'f (rn) , (3.20) 

iQ-ySj'l'k{rrn)^i^f{rn) = C (^n) ® (r^) ■ 
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where r^, > r„ (m > n) . 

By analogy with (3.6) we can define the "^-matrix as hnear map TZ : ^ (8) ^ — > ^ 

mnr) ® V^f (rO) ^ n^^J-i^^ f') ^lif) ® ^li^') = -^fptlir, r'|f, f') ^kif) ® <(r ), 

(3.21) 

where we imply summation over values ri, r^, r^, ...,rN of the discrete variables r,r\f,f' 
and X is the diagonal matrix 



lZij\r, r'|f, f') = X;,^,^;. <5(r, f) 5(r', f') 5(r, r') - - J^T^ S{r, f) 5{r\ f') (5(r, r') - 1) . (3.22) 
Let us require that the 7?.-matrix satisfies the quantum Yang-Baxter equation 

^12^23^12 = ^23^12^23, (3.23) 

and the Hecke relation 

-R? -{J - x)n - J = ( J + J • 7^) • ( J - i • 7^) = 0, (3.24) 



where J = J^p^Hr,r'\r,r') = J^j^p^,5{r,r) 6{r',r') is unit matrix, J' is diagonal matrix 



r a/3 k I 



-yptj 



jTi\r. r'lr, f) = 5(r, f) 5(r', f') 5(r, r') - I J.f^J 5(r, f) 5(r', f') (5(r, r') - 1) (3.25) 

and J = X ■ J . 

The 7^-matrix satisfying the Hecke equation (3.24) has the following block structure: 



Qf^f^ S{r, r') 6{r', r) 0"' + Q^^'^^; 6{r, r') 6{r', r) Q^'\ 



Kptlir, r'\r, f') = kf^^^ 5{r, r) 5{r', f') S{r, r') + (l - T^) J.f.f 5(r, r) S{r', f')e-'+ 

Xf^ Xf^' ^\ Oi^' r (3 26) 

where the f2-marix, and 6'"''' were defined in (3.10) and (3.11) respectively, S{r, f) is a discrete 
5-f unction. 

In the product ^ (g) in (3.21) we imply the following ordering: 



(3.27) 



With this ordering the T^-matrix has block structure and, for example, for two arbitrary sites 
r„ and (r„ < r^) the block "^(nm) of the 7^-matrix acting on the sum ^if^ ® ^'^^)©^'[^) ® 
^l'')©^!'') ® ^l'')©^!^) ® ^1^) has the form: 

' R " 

{l-l-^)Js Q 

^(nm) = Q Q-1 ■ {^.28) 

^ 
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By analogy with (3.14) and (3.15) from (3.24) we can define the quantum {l,q,s)- 
antisymmetrizer 

i=(J-X•7^), (3.29) 
and the quantum {I, q, s)-symmetrizer 



S={J + J-n) 

which are orthogonal projectors 

i.5 = <s-i = o, {Ay ■{i+j)-A, 



(3.30) 



{Sf = J^-it + J)-S. (3.31) 



The quantum symmetrizcr S permits to define the commutation relations for the lattice 
(/, q, s)-Grassmann bispinor field ip!^{r) {q^ 7^—1) 



5.(V'®V) = 0, or i^t{T)®i^^^{T') + {J-nr^^^^^{Ty\f,f')rk{f)®ri{^^^ (3.32) 

Let us require that these commutation relations are covariant with respect to the transfor- 
mation of lattice field il^f{r) by a quantum matrix A e SLq{2,C) and ip{r) is spinor with 
respect to SLs{2,C). These requirements lead to the following solution of the quantum 
Yang-Baxter equation (3.23) 



where Qg is defined in (3.17) and 



(3.33) 



Qs — -^s — 



s 







1 






1 



s 



Substituting this solution in (3.32), we obtain commutation relations 

€(rm) ® infirm) = -{Js ■ B)t'^ ^l{rm) ® <(r„) , 

C('^m) ® C(^n) = -lqs ^r(r„) ® (r„) , 
C('^m) ® (r„) = -/ s (r„) ® (r^) , 

i^tirm) ® ^,"(r„) = -Iqijjfirn) ® ^r(r„) , 



^^(rm) ® V'f (rn) = -Iqi^tirn) ® V'" (^m) 
where Tn < , i > j and a> (5. 



ls{q-q-')^r{rn)®i^Hrm), 
lq(s-s~^)'ijj^(rn)®'il^r(rm), 
/(s-s-i)^f(r„)®^f(r^), 

l{q-q-')ipUrn)®ipr{rm), 



(3.34) 



12 



4 Integral calculus for the lattice (/, g, s)-Grassmann 
field 

At first let us briefly recall the definiton of the G'Lg(2, C)-covariant generahzation of the 
Berezin integration over the quantum vector space ^^'^^ [24]. 

Let be space of the hnear functional on ^^^^ with dual nilpotent basis {J^'a} 
{{d'^ccY — 0) which maps ^^^^ in complex numbers: 

(^*(2)^^(2)) 

and 

(JVa, V'O = / dll^all^^ = 5i {dll^^, C)^ j dilJa C = 0, (4.1) 

where c is a complex number. 

Let us define the action of tensor product of the linear functionals by the following 
relation: 

V'^^V'") = = / d'^M'^-'^'^^V = 815P,{5^,-1). (4.2) 

Using this definition and the following ansats for commutation relations for dif^a'- 

di/jai <S> dip/s^ = X di/jp^ (g) JVai , 
where cti > we obtain (no summing) 

j di^oci (g) d^p^ij"^ (g) ^ 2 = -^y" d^Jai d^jp^jjj'^^ (8) -0"' = 
—qx j dijji3j^ (8) di/jaii/^^^ V'"^ = —qx. 



Prom here x — —- and 
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or in covariant form 



#„ ® + g (Rq) l^p d%l)^ (^diljp = 0, (4.3) 

where {Rg)ap = Rgja determines commutation relations for matrix elements of inverse 
matrix i"^ e GLg{2, C) 

iK) t (^~') I (^"') e = (^-') " {A-') I {K) II ■ (4-4) 

Hence, it follows: 

if V^f = ® V^'' , then d-^)^ = (^"^)^ ® d^fi- (4-5) 
Using (4.1) and (4.2), one can define the GL^ ( 2, C) -invariant Berezin integral on ^(^): 

J dipi ^dip2ip^ = 1- (4.6) 
13 



It is not hard to show that this integral is invariant with respect to transformations (4.5): 

JV^i (g) dip2 ijj^ ^ip"^ = detg{A~^) detgA J #i (g) #2 ® = 1 ■ (4.7) 
Using (4.2) and (4.6), one can calculate the integral 

#i#2^"^'^ = (4.8) 



where 



and for brevity the notation of tensor product is omitted. 

Now let us consider the integral calculus for the lattice (/, q, s)-Grassmann bispinor field 
ip = {ip°'{r)} which satisfy commutation relations ( |3.34| ). By analogy with (4.1) let us 

define a space \1> of the linear functional on \E' ( |3.19|) with dual nilpotent basis {dipl^{r)} 
([rf^^(r)p = 0): 

and 

(ciC(r),^f(r')) = / d^|Jl,{r)^|J^{r') = 6^J]6{r,r') {d^P\{T),c) = [ rf>;(r)c = 0, (4.10) 



where c is a complex number. 

Elements of the dual basis satisfy commutation relations 

#f (r) ® #f (r') + {J ■ n') ;;^f/(r, r'|f , f) d^PUf) ® d>f (f') = , (4.11) 



where T^'-matrix has blok structure ( ^.28[ ) and contains matrices (R^) = Rgp^ and 



( f?'\ kl _ jDlk 

Let us define the GLq{2, C)-covariant integration measure Vtfj for Berezin integral over 
the quantum vector space 

V^jj = d^/jlin) d^plin) d^lin) d^jlinWAr^) di^lir^) d^lir^) d^l{r2) ■ ■ ■ 

di:\{rr,) d^l{rN) diP^ir^) d^lir^,) . (4.12) 

By analogy with (4.8) and (4.9) we can define GLq(2, C)-invariant generalization of the 
Berezin integration for the lattice (/, q, s)-Grassmann field 

V^|J #^^(ri) #/(n) #/(ri) C/(ri) ... <^(r^) <^(r,v) <^(r,v) <^(r^) = 

01020304... Si52S-s54/ \ (AA'^\ 

Here the tensor e is defined by the following rules: 

1212... 1212/ N _ . 

1 122... II22I ' 1) 'NJ — -1- , 
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and the remaining components are defined by the coefficients appearing in the Lh.s of the 
relation 

- ^fA^N) <^(r^) <^(r^) <(r^) = 

r^) V^i(ri) ^l{r,) ^l{n) ... V^l(r^) ^l{r^) ^l{rr,) ^l{rr,) 

after reordering the l.h.s. to the r.h.s. using commutation relations (3.34) for each specific 
set of values of the indices. 

Now let us consider the one-site Gaussian integral over lattice (/, q, s)-Grassmann field. 
In order to calculate this integral, one uses the following definition of (g, s)-Pfaffian (which is 
an extension of the definition of usual Pfaffian proposed in [27] to g-deformed commutation 
relations). 

Consider the quadratic form 

w= all ^l + all ^2 + 2 ^2 + 2 ^? ^2 + 

a, 13, i, k 

where the matrix elements chosen to be commutative and the matrix a has the form 








"ll 


"l2 


"l2 









"l2 


22 
"l2 


a = 









"2 2 






sq "12 


~g "22 






Here the lower-triangle matrix elements are determined using relations ( p.l8| ). Let us define 
(g, s)-Pfaffian of d through 



Then, taking into account commutation relations (3.18), we get 



D J! ^ I 2 4\ 12 12 /I I 2\ 11 22 I 2 12 21 

Pfiq,s){a) = -{1 + q S )"ii"22 - 2^1 + ^ )"l2"l2 + 5? "I2"l2 

Using (4.13) for the one site integral, it is immediate to show that 

J dipl dtpl dipl di)l exp M ^ 'hk C [ = ^/(g,s)(")- 
Now let us show that, at g = ±1, s = ±1, 



Pfiq,s)ia) = ^Jdetib) 



(4.14) 



(4.15) 



(4.16) 
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where det is the usual (g = 1) determinant. To determine the matrix b we consider {q, s)- 
bispinor v, which satisfies commutation relations 

and require 

ivir = {vlf = {vlf = {vlf = 1 and 

In order to define these commutation relations we use quantum (g, s)-antisymmetrizer (3.15). 
Then, at g = ±1, s — ±1, i/^f = '4^t components of usual Grassmann bispinor 

Making use of the substitution ■^f — > ■^f = into integral (4.15), one gets 

vl v\ v\ v\ J di^l di^l di^l di^l exp 1 1 Kk # } = Pfi,,s){^) ■ 

a,l3,i,k ) 

where = (here summation is absent) . 

Calculating this integral, we obtain the connection between the usual Pfaffian of the 
matrix b and {q, s)-Pfaffian (4.14) 

vlvlvlvlPf{b)^Pf^,,s){a). 

This relation yields 

iPfi,,s){a)f = {Pf{b))'^det{b). 

Hence, we obtain (4.16). 

The {I, 5)-fermionization of the two-dimensional Ising model might be done using the 
lattice (/, g, s)-Grassmann field at q = I = —l,s — 1 (see the following section). In what 
follows, we just consider this case, although to compare with the case I — q — 1, one 
sometimes has not to fix the value of /, q. 

Note that for g = / = —l,s = 1 the (/, g)-Grassmann field '4'^{r) is real. To show this, 
we define the real field as V^f (r) — i^fir), where the bar denotes the Hermitian conjugation. 
Thus, we have 

(V'f(r)V'f(rO) = i^lir') ^(r) = V'f (r') C(0 ■ 
On the other hand, the Hermitian conjugation of fifth relation in (3.34) gives the relation 

The consistency condition for these relations leads to the restriction |/| =1. Using similarly 
arguments for (3.15), we get |g|^ = = 1. 

In further calculations we simplify the notation of the lattice (Z, Q')-Grassmann bispinor 
field 

V^(r) = {^\{r), ^l{r), ^l{r), ^^^(r)) = {^l, ^t) . 

and, as a corollary of (3.15) and (3.34) &X, I — q — — l,s = 1, its components satisfy the 
commutation relations 



•'fc 5 



0. 
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{^l , i^l^} = , = {C , C'} = , (4.18) 

where a = 1,2, 3, 4. 

Using commutation relations (4.17), (4.18) and (4.12)-(4.15), it is easy to calculate the 
one-site integral, which we need in the following section 

/ dl/jr exp I bik i^i^ i^r, + h i^l, i^rs ^t, } U expi^J^J = 
[i<k ) i=l 

= b4 + Pf{l,q)ibik) + bi2 + hi + &24 + ^13 - ^14 - h23, (4.19) 

where i , A; = 1, 2, 3, 4, / = g = — 1 and 

y JVr = j #lri #2r2 ^Srg #4r4 = , 
Pf{l,q){hik) = ^'12 &34 + &13 &24 - &14 &23- (4.20) 

5 (/, g)-Fermionization of the two-dimensional 
Ising model 

In this section we discuss how to express partition function of the two-dimensional Ising 
model through functional integral over the lattice (/, g)-fermion field {I = q = —1). The 
method allows one to derive the Lagrangian of the lattice (/, g)-fermion field theory describing 
the behaviour of the Ising model nearby the critical point. Our method works for the two- 
dimensional Ising model of general type with the only restriction of Z2-invariance of the 
plaquet statistical weight. In this case, the Hamiltonian /i((Tfl) (energy per plaquet) includes 
not only the pairwise nearest-neighbour interactions, but also four-spin interaction 

4 

h{(rR) = -J^-^ij^^- -^IK^ i,j = l, ..,4, (5.1) 

i<j i=l 

where (j^ = ±1 is the Ising spin, R = {x,y), x = 1, .. ,n, y = 1, .. ,m, denote the 
plaquet position on the lattice (as before we use index r for numbering lattice sites), Jij are 
the parameters of the pairwise interactions and J is the four-spin coupling constant 
(see Fig.) 

Note that for the Ising model with only the nearest-neighbours interactions we have 

Ji — J34 -I J2 — J24 ; J — J\2 — Ji3 — J24 — Ji4 — 0, 

where Ji and J2 are the vertical and horizontal coupling constants, respectively. 
The Ising model Hamiltonian on the whole lattice has the following form 

h^Y: hM = - E(E Jv< < + (5-2) 

R R i<j i=l 
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Y 

Here x and y denote the unit vectors along the horizontal X and vertical Y axes, respectively. 
This model corresponds to the eight- vertex model in external field [28], whose path integral 
obtained at the end of the section is not calculable due to the four-fermion interaction. 

Partition function for the model with Hamiltonian (5.2) can be represented in the form 
of the product of plaquet statistical weights 

where cr^ is the Ising spin on site r and one gets at i? = r 

<^fl((7r , CTr+y \ (^r+x , f^r+x+y) = ^^('^fl ^'^II'^r, '^r) = ■ (5-4) 

It is easy to show that statistical weight (5.4) can be presented in the following form 

4 

uj^{a) = p(l + ^ aik(Ti (tI + ai [] 4)> (5-5) 

i,k i=l 

where 

Oiik^< crl^a'^ >h, q;4 =< crjj (7^ C7| (7^ >fe , p = 2"*^ ^ (^r{(t) , (5.6) 

are expressed through the initial coupling constants Jik , J and < /(cr) >h denotes averaging 
over the plaquet statistical weight 

< /(<7) >h= E /(<7)^h((7)/ E 

o-i,.,o-4 0-1,., 0-4 

Note that, choosing the statistical weight in the form 

4 

ujR{a) = exp(/Ci al ^ + /C2 4 + /C <) , 

i=l 
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where /Ci = /3Ji , /C2 = /5^2 , JC = PJ, one gets the eight- vertex model in zero external field, 
and, for the Ising model with only the nearest-neighbours interactions, one obtains 



ujR^a) = exp(/Ci cr^ al + IC2 (^l o"^). 

In the last case, using (5.6), it is easy to calculate the coefficients in (5.5) 

p = cosh /Ci ■ cosh /C2 5 ^34 = h , 0^24 = ^1 , 
a23 =tit2, a;i4 = 012 = 04 = 0, 



(5.7) 



(5.^ 



where ti = tanh/Ci and ^2 = tanh/C2. 

To obtain partition function (5.3) as a functional integral over the lattice {I, g)-fermion 
field, the key points are: commutation relations (4.17) and (4.18), definition of integration 
rule (4.13) over (/, g, s)-Grassmann field at / = g = — l,s = 1 and representation (5.5) for 
the plaquet statistical weight. The latter one can be written using (4.19), (4.20) in the form 
of integral 



/I 2 I 3 4 \ 



= p Jdipu exp [ J2i<k aik cri cr^ -^i 
«4 crl o\ a\ a\ ^\ ^\ ] nti exp [ 



(5.9) 



where (h\}j^ = dtpi ^ dtp2 r d^3 Rd^4R- 

The parameters in (5.5) and (5.9) are expressed through each other as 

ai2 = ai2, = 013, 0:24 = ^24) Ct34 = O34; Ct24 = O24, 

ai4 = -ai4, ^23 = -^23, ^4 = 0-4 + Pf-i{aik), 



and (/, g)-fermion field components ip satisfy commutation relations ( [4. 171) and ( |4.18|) . 

In a similar way, one can write the representation of the statistical weight product in 
partition function (5.3) 



R •' 



(5.10) 



where 



Cr {(Tip) = ar (Tr+y{ai2 Ipl + 034 ^J^+y) + ^^r C^r+£(ai3 V^r V'r+x + «24 i^r i^t+x) + 

+ar ar+x+y ^14 ^t+x+y + ^r+£ (^r+y ^23 ^^r+y ^r+x + 



+ 04(0",. ar+x (Tr+y (Jr+x+j/) (V'r "^^r+y ^r+x ^r+x+y)^ 

pv^=n(^^>r-). 



(5.11) 



Since not all of the multipliers in Vr and Or are commutative, it is necessary to indicate 
their specific ordering. In (5.10) we imply the following ordering 



n rn 
r y=l x=l 



(A,l -^2,1 • ■■'Pm,l)X 
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{Vl2 ■ V2,2 ■ ■ ■ Vm,2) ' ' ' (7^1,n ' 7^2,n ' ' ' Vm,n) ■ (5.12) 

Relation (5.10) is proved by termwise integration over the field components ipl , , 

r+x+y ■ 

Product (5.12) can be re-odered in the following way 



where 



Qr = e^' e^' e^" e^" . 

Let us emphasize that, in the course of this rc-odcring, we do not rearrange the non- 
commuting exponentials, for example, exp(-(/'^) and exp(-(/'^') or exp(-(/'^) and exp(^^,). One 
can easily check it, writing down Vr in (5.12) explicitly through the corresponding exponen- 
tials in (5.11). Right in order to do this re-odering, we need exotic commutation relations 
(4.17), (4.18) for the lattice (Z, q')-fermion field I — q — —1. Using this commutation 
relations and the nilpotence of the components ^0^, one transforms Qr to the form 



Qr =(l + ^V';)exp[£, 



(0) 
r 



i=l 



where 



As a result, one gets for integral (5.10) 



R 



exp 



n (i + EV';)exp[4°) w 



(5.13) 



The next obvious step is to replace the variables 

which leads to vanishing the terms Uj. cr^' in integral (5.10). Therefore, the dependence of 
the Ising spin is localized in the multiplier (1 + (jRYli=i Vl) 

i=l i=l 

Cr{a^P) = Criv), Cf{^)=Cf{^). 

Now the summation over ar in partition function (5.3) is independently done at each lattice 
site 

(Tr=±l 1=1 

and the "bad" multiplier in product (5.13) vanishes. Since terms in {(p) and Cr (9?) are 
quadratic over the fields (fr (with an additional quartic term provided 04 7^ 0) and, therefore, 
commutative, one obtains 



Dip exp 



X:(4°^ {^)+Cr{^)) 



(5.14) 
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where 

+ai4 (pi (pt+x+y + «23 pI Pr+x+y + «'4 pl pl+y Pr+y vt+x+y ■ 

At a4 — 0, one can represent functional integral (5.14) as Gaussian integral 

1 



Vp exp 



where 



D 



1 





ai2 V_ 



1 + ai2 




ai3 + 023 + Va 

ai4 V_y + a24 + V-a; 



ai3 - V-^r 
023 V_y + V_ 



1 + V_„ 



(5.15) 



ai4 Vy + V_ 

024 ~ 

1 + a34 Vy 





This integral is expresed through (/, g)-Pfaffian of the matrix D. Similarly to the case of 
usual fermionization of the two-dimensional Ising model (1.2)-(1.4), this matrix is linear in 
the one-step shift operators Va; , Vy. In fact the accounting of the boundary condition in 
(5.15) gives the (Z, g) -deformed counterparts {I — q — —1) of (1.2), (1.3). For the sake of 
simplicity, we consider only Z^^, where index "A" again denotes the antiperiodic boundary 
conditions along the axes X and Y. When turning to the momentum representation, the 
momentum components pi,P2 take half- integer values. 

In what follows, one restricts himself to considering the two-dimensional Ising model with 
only the nearest-neighbours interactions. In this case, the coefficients aij are written in (5.8) 
and functional integral (5.15) in momentum representation can be rewritten in the following 
form 

Z={2pY^ j P(^exp{^ v\-Pi.-P2)D^KpuP2)pKpi.P2)}. 



(5.16) 



where 




1 



1 



e~'P' - tit2e 



-e -'^^ e 
tit2e-'P^ ti-e 



ipi 



IP2 



ti + e 



-ipi 







1 + ^26- 



1 + t2e'P^ 



and 



Vp= Yi dRepi{pi ,p2) dRep2{pi ,^2) dRepsipi ,^2) dRepi{pi ,^2) 

P<>0 

dlm(pi(pi ,p2) dlm(f2ipi ,^2) dlrrnpsipi ,^2) dlmp^^ipi ,^2) ■ 

At q — —1, s = 1, using (4.16) for the calculation of the {q, s)-deformed Gaussian functional 
integral, one immediately gets integral (5.16). As a result, one obtains known expression 
(1.6) for Z^^. 
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Thus, (5.16) connects the partition function of the two-dimensional Ising model with 
functional integral in the lattice free (Z, q)-ieicimon field theory (I — q — —1). The coincidence 
of the results of calculations (5.16) and (1.8) allows one to establish the validity of the 
definiton of the Berezin integral for the lattice (/, g, s)-Grassmann field, proposed in (4.13). 

Note that in paper [29] it was demonstrated the connection between paragrassmann 
algebras and the some representation spaces of quantum matrix group I/q(2, C) with q being 
a (p + l)-root of unity where p is integer (^*'+^ = for paragrasmann 6). In our case, using 
definition of the grassmann anyonic field given in [30 ], the commutation relations (4.17)- 
(4.18) for the lattice (/, g, s)-fermion field {I = q = —1, s = 1) can be considered as definition 
of the lattice two-component paragrassmann anyonic field with / = e*'^^, {y = 1) and p = 3. 
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